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Abstract 

We compute the asymptotics of a block Toeplitz determinant which arises in the 
classical dimer model for the triangular lattice when considering the monomer-monomer 
correlation function. The model depends on a parameter interpolating between the 
square lattice (t = 0) and the triangular lattice (t = 1), and we obtain the asymptotics 
for < t < 1. For < t < 1 we apply the Szego Limit Theorem for block Toeplitz 
determinants. The main difficulty is to evaluate the constant term in the asymptotics, 
which is generally given only in a rather abstract form. 

1 Introduction 

This paper is devoted to computing the asymptotic behavior of a certain determinant that 
arises in the study of the classical dimer model on the triangular lattice. In jB], Fendley, 
Moessner, and Sondhi study the monomer- monomer correlation function p( mrn \n), which 
can be expressed by a determinant of a block matrix 

P {mm \n) = l -^M~ n , M n = 

where TZ is an n x n matrix with entries depending on the difference of their indices, and Q 
is an n x n matrix with entries depending on the sum of the indices. The entries also depend 
on a parameter t which interpolates between the square lattice (t = 0) and the triangular 
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lattice (t = 1). The authors computed this correlation function numerically for t = 1 as 
the size of the matrix increases and found that it converged to a constant value of around 
0.1494 . . .. For all the details of the dimer model computation, the reader is referred to j^]. 

The main purpose of our paper is to compute the determinant asymptotically for all 
values of t between zero and one. In fact, we will be able to determine the asymptotics for 
all complex parameters t with Re (t) > 0. This covers all physically interesting cases except 
t = 0. 

Our method is to convert the determinant of the block matrix M n that arises in the dimer 
problem into a determinant of a block Toeplitz matrix and then to find a way to explicitly 
compute the asymptotics of the block Toeplitz determinant. 

To be more specific the dimer matrix ((TJ has n x n matrix blocks 1Z and Q whose entries 
are given by 

TZ jh = 2(-l)^-^R k _ 3+1 + 9(j-k^- k - 1 
Q jk = 2i(-l)^ +k ^Q n+1 _^ k , 

(1 < j, k < n) and where the expressions R k , Qk, and 6(k) are defined by the following. 
For even k 



and Qk = 0. 
For odd k 

and 



IFF cosycos(kx + y) 

k = sT~2 / / 2~^ 2 T - 2 27 , \ d y dX 

6ir z J_ n y _ 7r cos^ x + cos z y + t z cos z {x + y) 



^ rn r-K COS X COS kx , , 

Q k = Q~^2 / / 2 7T2 27 -^ d V dx > 

57r z J_ w J_ T cos z x + cos z y + V cos z {x + y) 
t F F cos(x + y) cos(kx + y) , , 

k = \ \ 2 i 2 rT2 27 , \ d V dx - 

&7T Z J _ n J _ 7r cos^ x + cos^ y + t z cos z {x + y) 



The expression 6(k) equals 1 for k > and otherwise. 

In the Appendix we show that the determinant of the dimer matrix M n is the same as 
the determinant of an n x n block Toeplitz matrix. A finite block Toeplitz matrix is one of 
the form 

T n {4>) = (<f>j-k), < j, k < n - 1, 

where 

1 f 27T 

' (e ix )e~ ikx dx 



k 2tx 



o 



are the (matrix) Fourier coefficients of an iV x iV matrix valued function defined on the 
unit circle T = { z G C : \z\ = 1 }. In the dimer case we have N = 2, and for < t < 1 we 
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obtain 

(c d 
I , ( 2 ) 
d c 

where 

^ = (tcosx + sin 2 x) ^ 

(t — e~ lx )Vt 2 + sin 2 x + sin 4 x 

die") = . 4 = ■ (4) 

VI + sin x + sin x 

Here and in what follows d{e lx ) = a(e~ lx ). 

In the case < t < 1, we are able to find the asymptotic expansion of the determinant of 
T n ((j)) by using the Strong Szego Limit Theorem. The main difficulty is that the constant term 
in the asymptotics for block Toeplitz matrices is generally only given in a rather abstract 
form. Thus we are required to find a more explicit expression for the constant, which 
fortunately, in this case, can be done. For this purpose we will use an expression for the 
constant which holds in very special cases and which is due to Widom We will also use 
another determinantal identity which is perhaps new and by the help of which more general 
(yet special) cases can be reduced to the situation covered by Widom. We will also give a 
new proof of Widom's result using the Borodin-Okounkov-Case-Geronimo identity PQI21IZj- 
All these results are described in the Section El and are of independent interest. 

The application to the dimer case still requires a lot of elementary, yet tedious compu- 
tations which will be done in Sections 0] and 03 The final result is that the limit of the 
monomer-monomer correlation function 

P (mm) (oo) := lim P^ mm \n) (5) 

n— >oo 

exists and is given by the formula 



p(™» \oo) = - J . (6) 

At this point we have proved this in the case < t < 1. 

In the case t = 1, we will also obtain this asymptotic formula, which agrees with the 
numerical calculation In this case we cannot use the Szego Limit Theorem because the 
function c has a singularity. This singularity is due to the second term in the definition of 
the entries of 1Z in terms of Rk- The entries of M n are well-defined and analytic for t G C 
with Re (t) > 0. For these values we will actually be able to prove the same asymptotic 
formula as above. This will be done in Section El The idea is to transform the matrix T n ((f>) 
into a matrix which is a finite block Toeplitz matrix plus a certain perturbation (which is 
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of fixed finite rank). A generalization of the block version of the Sego Limit Theorem then 
gives the expected asymptotics. The statement and the proof of this generalized version is 
given in Section EJ and is also of independent interest. 



2 Determinant identities 

We begin with some preliminary facts about Toeplitz operators and Toeplitz matrices. Let 
<f) g L°°(T) NxN be an essentially bounded N x N matrix valued function defined on the unit 
circle with Fourier coefficients (f>k G C NxN . The Toeplitz and Hankel operators are defined 
on £ 2 (Z + ), Z + = {0, 1, . . .}, by means of the semi-infinite infinite matrices 

= 0<j,k<oo, 
H{<t>) = Oj+fc+i), 0<j, A; < oo. 
For 0, G L°°(T) NxN the identities 

r(#) = T(<j>)T(iJ>) + H{<f>)H$) (7) 
H(W>) = T{<j>)HW)+H(<j>)T$) (8) 

are well-known. It follows from these identities that if ip- and ip + have the property that all 
their Fourier coefficients vanish for k > and k < 0, respectively, then 

T{^_U + ) = T(^_)T(0)T(^ + ), (9) 
H(^ + ) = T^)H(<P)T^ + ). (10) 

In the following sections we compute an explicit expression for one of the constants that 
appears in the version of the Strong Szego Limit Theorem for matrix-valued symbols, for 
certain symbols having a special form. In order to state the classic Strong Szego Limit 
Theorem for matrix- valued symbols, let B stand for the set of all function G L X (T) such 
that the Fourier coefficients satisfy 

1/2 



E w + ( E i*i -i^i 2 ) <o °- ( n ) 



With the norm (|TT|) and pointswise defined algebraic operations on T, the set B becomes a 
Banach algebra of continuous functions on the unit circle (see, e.g., [21 EI])- 

Theorem 2.1 (Szego- Widom [10j) Let <fi G B NxN and assume that the function det 
does not vanish on T and has winding number zero. Then 

det T n (<t>) ~ G(0) n £(0), as n -> oo, 
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where 

G{<f>) = exp Q- log det cj){e ix )dx} 

and 

E(<f>) = detT{(j>)T{(f)- 1 ). 

In connection with this theorem, let us explain why the definitions of the constants G(<f>) 
and E((f)) make sense. First of all, one can show that each nonvanishing function in B with 
winding number zero possesses a logarithm which belongs to B. Hence G(4>) is well defined 
as the geometric mean of det <fi. The constant E((f>) has to be understood as an operator 
determinant. In fact, we have T(</>)T(</> _1 ) = I — H^cfyH^^ 1 ) where the product of the 
Hankel operators is a trace class operator. Observe that we have defined B is such a way 
that Hankel operators with symbols in B are Hilbert-Schmidt. For general information about 
trace class operators and operator determinants see, e.g., jHj. 

In the scalar case (N = 1), there exists an explicit expression for E(<p) given by 

oo 

E{4>) =exp ( ^ A; [log 0] fc [log J. 
fc=i 

In general (N > 2) a more explicit expression for E((f>) does not exist apart from very special 
cases. For more information about related results we refer the reader to [2|. 

One of the few results concerning an explicit expression for E((p) in the matrix case goes 
back to Widom who showed that if <pk vanish for k < n or for k > n (for some fixed n), 
then E(<f>) = det T n (0 _1 )G(0) n (see 0, Theorem 5.1). Thus the constant E((p) is reduced to 
the computation of a finite determinant, which for small values of n is very computable. If 
one can somehow reduce the computation of our E(<f>) to this case, then an explicit formula 
for E((j)) is possible. The next theorem and its corollary will facilitate such a reduction and 
since we cannot expect to do it except for very special cases, it requires special assumptions 
on the symbol. 

Proposition 2.2 Let R G B NxN . Then the operator determinant det(T(e R )e~ T ^) is well 
defined. Moreover, /(A) = det(T(e Ai? )e _T ^ A ^) is an entire function. 

Proof. This is proved in Sect. 7]. □ 
In what follows let In denote the N x N identity matrix. 

Theorem 2.3 Let Q G B NxN , trace Q = 0, and a G B. Then 

det(T(e alN+Q )e- T{alN+Q) ) = exp (— trace #(a)#(5)) ■ det(T(e Q ) e - r(Q) ). 
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Proof. This proof is modeled on the ones given in 0J|3] and a more detailed account of why 
the various determinants and derivatives are defined can be found there. We give a only a 
sketch here. Let R = al^ + Q and define 

/(A) = det(T(e Ai? )e- AT ^). 

The function /(A) is analytic on C and it is nonzero if and only if the Toeplitz operator 
T(e XR ) is invertible. In particular, /(A) is nonzero except on a discrete subset of C. Using 
the fact that (log det F(X))' = trace F^A)^ 1 (A) we have that 

/'(A)//(A) = trace (t (Re^^T -1 (e XR ) — T(e XR )T(R)T^ 1 (e XR )^j . 

This implies 

/'(A)//(A) = trace (r-\e XR )T{Re XR ) - T(i?)). 
Differentiating again we have that 

(/' (A)//(A))' = trace (r-\e XR )T(R 2 e XR ) - T- l {e XR )T{Re XR )T-\e XR )T(Re XR )^j . 

At the points where /(A) 7^ 0, the symbol tp x = e XR has the property that it factors into a 
product 

?P X = tp x ip x 

such that the factors ip- and ip + as well as their inverses belong to B NxN and have Fourier 
coefficients that vanish for k > and k < respectively. Thus, by Q, 

T(e XR ) = T(^ A )T(^ A ). 

This yields 

T- 1 (e A/? )=T(^ A )- 1 )T((^ A )- 1 ). 

After simplifying we obtain 

(/'(A)//(A))' = trace (t((^ a )" 1 j R^ a ) - T 2 ((^)- 1 ^)) = tx a ce(H(g)H(g)) 

where 

q = (^ x )- l R^ x = aI N + (^ A )~ 1 Q^ A . 

Now define 

g(X) = det(T(e AQ )e- Ar(Q) ) 

and compute the above expression once again with R replaced by Q (i.e., put a = 0). Then 
the difference 

(/'(A)//(A))' - Q/(A)/<7(A))' = N ■ trace H(a)H (a). 
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The reason for this is that if we expand 

trace H(g)H(g) 

we have four terms. A term of the form iif (aIjv)-ff(('0-)~ 1 Q'0-) has trace zero since each 
matrix Fourier coefficient of (ip ) l)~ 1 Qip > l_ has this property and multiplication by a matrix 
coefficient of al^ corresponds to a scalar multiplication. The term 

trace 

cancels with the corresponding term for g since the corresponding functions if>* (obtained 
from e AaJ iv+ A( 2 an d e A< ^, resp., by factorization) only differ by a scalar function. Finally the 
"iV" comes from the fact that aijv is a N x N matrix. 
Thus 

log /(A) -log g(X) 
has a constant second derivative and satisfies 

(tog/-logp)(0) = and (f/f)(0)-tf/g){0) = 0. 
From this the theorem follows easily. □ 
Corollary 2.4 Let Q G B NxN , trace Q = 0, and a G B. Then 

E(e alN+Q ) = exp (n ■ trace H(a)H(a)J ■ E(e Q ). 

Proof. We can write 



E{e alN+Q ) = detT{e alN+Q )T(e- alN - Q ) 

= det T{e alN+Q )e~ T{alN+Q) ■ det e T{alN+Q) T(e 



air, 



and apply the previous theorem twice. □ 

The following theorem and proposition yields an alternative proof to Widom's result in 
the case where either the positive or negative Fourier coefficients vanish for all but a finite 
number of indices. We give a proof that is based on the Borodin-Okounkov-Case-Geronimo 
(BOCG) identity PQE]. Widom's result is stated in our Theorem 12.71 

Proposition 2.5 Suppose that if) G B NxN such that both T(if>) and T{if>) are invertible on 
(P(Z+)) N . Then 

detW 1 ) = ^^^(l-H(z^)T-\^)H(^)T-^)). (12) 
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Proof. The invertibility of T(ip) and T(ip) implies the invertibility of ip and the existence 
of a left and a right canonical factorization (in B Ny - N ) 

Proceeding as in the proof of the BOCG-identity one obtains 

detX^- 1 ) = G(%l)~ l ) n ■ detP„T(^ 1 )T(n_)T( M+ )T(i;: 1 )P ri , 

where P n is the projection {xk}^L Q i— > {a^o, • • • , x n -i, 0, 0, . . .} acting on £ 2 (Z + ). Using the 
identity 

det P n AP n = det(Q„ + P n AP n ) = det(Q„ + AP n ) 

= det A- det(A^ 1 g n + P n ) = det A ■ c\et(Q n A~ l Q n ) (13) 

in which Q n — I — P n and A is an invertible operator of the form identity plus trace class, 
it follows that 

det P n T(v+ l )T(u-)T(u + )T(vZ l )P n 

= det r(t;; 1 )r( M _)r( M+ )T(t;z 1 ) • det Q n T(v.)T{u+ x )T{uz x )T(v+)Q n . 

The first term on the right is easily seen to be det T^T^ -1 ) while the second term is equal 
to 

de\ 1 T{ Z - n v^uZ 1 )T{uZ 1 v + z n ) = det (J - H{z- n v_uf)H{u-_ l v +z - n )) 

where we use z~ n V-.uZ\ uZ v + z n = In. Now we substitute uZ 1 = ipu- and uZ 1 = u + ip and 
pull out the remaining factors from the Hankel operators (see p()jl ). which give the inverses 
of the Toeplitz operators. □ 

Proposition 2.6 Let ip £ qNxn ^ ^nfertt^/e. Then detT n (^ _1 ) = if and only if the 
operator 

( T(f) H^z~ n ) \ 
has a non-trivial kernel (or, equivalently, a non-trivial cokernel). 

Proof. The proof relies on the general fact that the kernel (resp. cokernel) of an operator 
PAP + Q is trivial if and only if the kernel (resp. cokernel) of the operator QA~ l Q + P 
is trivial, where A is assumed to be an invertible operator, P is a projection operator and 
Q = I — P. Indeed, this follows from the fact that both operators are equivalent, i.e., there 
exist invertible Si and S 2 such that PAP + Q = Si(QA^ 1 Q + P)S 2 - (One can also say 
that PAP and QA~ X Q are equivalent after extension. See (|13|) for the underlying line of 
argumentation.) 

We apply this statement in the setting A = M(ip~ l ) being the Laurent operator acting 
on (£ 2 (Z)) N , and P being the projection {x k } keZ ^ {xkjkZo on {£ 2 {%)) N ■ D 
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Theorem 2.7 (Widom [{J|) Let if) G B NxN be such that the function detip does not vanish 
on T and has winding number zero. Assume that ip^ — for all k > n or that = for 
all k > n. Then 

E^)=G(^) n detT n (?p- 1 ). (14) 

Note the above result holds also for n = when stipulating detTo(^ _1 ) = 1. 

Proof. The winding number condition implies that both T(ip) and T{ip~ x ) are Fredholm 
operators with index zero. Since one can show that T(-^ _1 ) and T(ip) are equivalent after 
extension, T(ip) is also a Fredholm operator with index zero. 

Hence if E(ip) ^ 0, then both T(ip) and T(^ _1 ) (and also T(if))) are invertible. Now 
Proposition I2.5I settles the assertion. 

If E{iji) = then (since T(^)T(-?/'~ 1 ) = /—compact), the product T^T^" 1 ) has a 
non-trivial kernel and cokernel. This implies that T(ip) has non-trivial kernel and cokernel, 
or, that T(?/> _1 ) (hence T(ijj)) has non-trivial kernel and cokernel. 

Consider the case that ipk — for all k > n. (The other case can be treated analogously.) 
Then the operator considered in Proposition 12.61 takes the form 

( T$) H(4,r n ) \ 

V o r(vo 

In the case where T(i()) has a non-trivial kernel, the operator above has a non-trivial kernel. 
However, if T(ip) has a non-trivial kernel (hence a non-trivial cokernel), the above operator 
has a non-trivial cokernel, too. □ 

3 A generalized Szego- Widom limit theorem 

When dealing with the case of parameters t ^ (0, 1) in Section |Ul we are going to employ a 
generalization of the matrix version of the Szego- Widom limit theorem. The reader who is 
not interested in the details of this case can skip this section. 

In the analysis of this case we are led to compute the asymptotics of the determinants of 
block Toeplitz matrices which are perturbed in some way by trace class operators. Although 
this generalization is very much straightforward, we give the details of the proof because 
our matrices depend analytically on the parameter t and we want to show that the limiting 
constant also depends analytically on t. This will enable us to identify the constant. 

The kind of sequences for which we are going to establish the limit theorem are the 
N x N block versions of the sequences described as follows. Recalling first the definition of 
the Banach algebra B of smooth functions given in the previous section, let T stand for the 
set of all sequences (An)^^ of n x n matrices which are of the from 

A n = T n (a) + P n KP n + W n LW n + C n 
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where a G B, the operators K and L are trace class operators on £ 2 , and C n are nxn matrices 
tending to zero in the trace norm. The set of such sequences {C n ) c £ =1 will be denoted by N ' . 
The operators P„ and W n acting on £ 2 are defined by 



P n ■ {x k }T=o {%o,Xi ■ ■ ■ ,x n _i, 0,0, • • •}, 
W n : {x k }^ =0 ^ {x n -!,x n -2, ■ ■ ■ ,x ,0,0, . . .}, 

and P n KP n and W n LW n are identified with nxn matrices in the natural way. 

Now we are precisely in the setting considered in j^j. (There is a change of notation, 
namely, that our B corresponds to S of The set T is a Banach algebra with algebraic 
operations defined elementwise and a norm 

||(A„)||^= ||a|| B + ||iqi + ||L||i + sup||C n ||i 

n>l 

where || ■ ||i refers to the trace norm. The subset Af is a closed two-sided ideal of T . 

Let gB NxN stand for the group of all invertible elements in the Banach algebra B NxN , 
and denote by QiB NxN the connected component of QB NxN containing the identity element. 
We remark that for a G Q 1 B NxN the constant 

( 1 f 2n 

G(a) = exp I — / log det a{e tx ) dx 

is well-defined Sec. 6]. 

In what follows we are going to consider elements (A ntt )^ =1 G J rNxN which depend 
analytically on a parameter t G fl. By this we mean that the map t i— > (A n ,t)5£Li £ j^ NxN is 
an analytic j^ NxN - valued function. If 

(Am)"=i = (r„(ot) + P„^P„ + W n Wn + C n>t )^ =1 G ^ xAf (15) 

this is equivalent to requiring that the maps t i— ► at, t i— > i^, t i— ► L t , and £ i— > (C nit ) are 
analytic (because of the definition of the norm in J rNxN ). 

Theorem 3.1 Let Q be an open subset of C. For each t G O let (A njt )™ =1 G J rNxN ; and 
assume that the map t G i— > (v4 njt )^ =1 G J rNxN is analytic. Moreover, suppose that 
a>t £ Q\B NxN where at is given by (li.5)). T/ien /or eac/i t 6 O £/ie /zrm£ 

detA n ,t 

Ei = lm (16) 



n^oo 



exists, the convergence is locally uniform on Q, and E t depends analytically on t. 
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Proof. By the Vitali-Porter Theorem on induced convergence |3J Chap. 9], it suffices to show 
that the sequence G(a t )~ n ■ det A Utt converges pointwise and is locally uniformly bounded on 

a 

Because of the assumption a t G Q 1 B NxN , we have a pointwise representation of the form 



a t = e b ^--e b « 1 b u . . . , b R G B N * N 

with R possibly depending on t. A moment's thought reveals that for each point t £ Q one 
can find a representation of the above kind on a sufficiently small neighborhood fi(io) of to 
such that bi, . . . , &r depend analytically on t. In what follows we will consider any fixed, but 
arbitrary point t <E Q and a corresponding O(to) ■ 
For t G fi(£o) consider 

B n ,t = e- T " (6l) • ■ ■ e~ T ^ (T n (a t ) + P n K t P n + W n L t W n + C n , t ) . 

The determinant of B Utt equals G(a t )~ n ■ det A n ^ since the trace of T n {b\) + . . . + T n (6^) 
equals n times the zero-th Fourier coefficient of trace(6i + . . . + &r) = log det a t . Moreover, 
since {B n ^) is built from elements in J^ NxN depending analytically on t, it also belongs to 
jtNxn an( j d e p enc i s analytically on t. 

Proceeding as in Prop. 9.2], one can show that the sequence B n j is of the form 

B n . t = P n + P n K t P n + W n L t W n + 6 n , t 

with K t , L t being trace class, (C ni t)™ =1 belonging to J\f NxN \ and with analytic dependence. 
Applying the Lemma 9.3 of [5 , it follows that 

lim det B nt = det (I + K t ) ■ det(J + L t ) 

pointwise. Since | det(J + A)\ < exp \\A\\i for trace class operators A, it follows easily that 
det B n t is locally uniformly bounded. □ 

One might think of giving another simpler proof of the previous theorem by considering 
the limit of the quotient 

<lrt l " ' - det(P n + T-\a t )P n KP n + W n T~ l (a t )LW n ) . 



det TJa t ) 



In order to make this work one needs the assumption that T n (a t ) is stable (i.e. uniformly 
invertible as n — * oo). It is known that this equivalent to the operators T(at) and T(at) both 
being invertible. Unfortunately, the assumption on a t made in the theorem only guarantees 
that these operators are Fredholm with index zero. 
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4 The basic computations for < t < 1 



The first step in the calculation of the asymptotics of the correlation function in the dimer 
model is to show how the block matrix considered by the authors in jH] has a determinant 
that can be computed using the theorems of the previous section. As mentioned in the intro- 
duction, in [0] the authors show the correlation is one-half the square root of the determinant 
of the matrix M n given in (|TJ. 

In the Appendix A we show how to first convert the above coefficients as Fourier coef- 
ficients of certain functions in the case < t < 1. From this we are able to identify the 
symbol of the block Toeplitz matrix with N = 2, and obtain that 



detM n = detT n (0) with 0= (17) 




where 



t cos x + sin 2 x 



<e lx ) = : ======, (18) 

[t — e~ lx ) v t 2 + sin x + sin x 

d ^ lx ) = m -t x . 4 • ( 19 ) 

Vt + sin x + sin x 

Recall that a(e ld ) = a(e~ ld ). Since we now have a block Toeplitz matrix our next step is to 
apply the Szego Limit Theorem. For < t < 1 the symbol is a smooth function. If we 
compute the determinant of the symbol we find that after simplifying, 

det <f) = cc-dd= (20) 

t 2 — 2t cosa; + 1 

and thus det <fi does not vanish and has winding number zero. In addition the geometric 
mean G(4>) is one since 

\ogdet(f)(e ix )dx = -— log((l - te ix )(l - te~ ix ))dx = 0. 



v ' 2tt„ ., 

Hence the Szego- Widom Limit Theorem implies that the asymptotics are given by 

det T n (<f>) ~ E(<f>), n -> oo, (21) 
and the limit of the monomer-monomer correlation function (jSJ) exists and equals 

p(™»)(oo) = iv^)- ( 22 ) 
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Note this agrees with the result from 6J that the determinants should approach a con- 
stant. Now we wish to apply Corollary 12.41 and Theorem 12.71 to identify the constant E(<f>). 
We write 

' V q 



a 



q p 



where 

p(e ix ) = (tcosx + sin 2 x)(t-e ix ), q(e ix ) = sinx(l - 2t cosx + t 2 ) } (23) 

and 

a(e ix ) = [f + sin 2 x + sin 4 x)~ 1/2 (l - 2t cosx + t 2 )-\ (24) 

The following lemma is needed in order to verify the assumptions of Corollary 12.41 

Lemma 4.1 For < t < 1, the matrix function <ft can be written as <ft — e al2+b ® with 
a,b G B, Q G B 2x2 , trace Q = 0. In particular, 

1 2 9 

a = — log(l — 2t cos x + £ ) + iri, Q = 

V ? ¥ 

Proof. Clearly, a and Q satisfy the stated conditions. Let us mention that b is given by 



and let us first show that b G B (if the logarithm is chosen appropriately). Notice that 
p + p 



(t cos x + sin x)(t — cosx), A = zsinxy (1 — 2t cosx + t 2 ) 2 + (t cosx + sin x) 2 . 

Thus for a(x) = — — A, we have lma(x) < for < x < ir, Ima(x) > for —n < 
x < 0, and both a(0) = t(l — t) and a(7r) = t(l + t) are real positive. Hence a(x) is 
a continuous function with winding number zero, which possesses a continuous logarithm 
which we normalize in such a way that the logarithm is real positive for x = and x = it. 
With this normalization the numerator in the expression for b becomes zero at x = and 
x = 7r, which cancels with the zero in the denominator A (as one can easily show). Hence 
beB. 

From the definition of b we further conclude that 

J* = e-V (E+E + A 
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We claim that 



This can be seen most easily by considering the product and noting that 

(P+^! _ A 2 = (t 2 + sin 2 x + sin4;r)(1 _ 2tcQSX + f 2) = e 2 V -2 

It follows that 

cosh(A6) = e-V ^t^, sinh(A6) = e"VA. 
It can be verified straightforwardly that 

= o- ^^/ 2 + , Q 2 = _(det Q) J 2 = A 2 / 2 . 

Hence 

ho i /, a x t sinh(6A) 
e feQ = cosh(6A)/ 2 + x — ( 2- 



A 

Combining all this yields 

e bQ = e-V ( P -±?I 2 + Q 



and thus <p = e ah+bQ . □ 

In the following computation the factorization of the function t 2 + sin 2 x + sin 4 x, which 
appears in the function a, will play an important role. We make a substituion z = e tx and 
can rewrite 

t 2 + sin 2 x + sin 4 x = -^{z A - 8z 2 + (14 + 16t 2 ) - 8z~ 2 + z' A ). 

Making the obvious substitution z = y 2 we can now factor 
y 2 - 8y + (14 + 16t 2 ) - 8y _1 + y~ 2 

= erva-^a-wa-ei!/" 1 )^-^" 1 ), (25) 

where £1 and £ 2 are defined by 

6 := 2 + /1 - 2 a/1 - t 2 + //, 6 := 2 - /i - 2 a/1 - t 2 - Ji (26) 

with 

/i := Vl - 4t 2 . (27) 
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Notice that |£i| < 1 and |£ 2 | < 1- Their inverses are given by 

= 2 + fx + 2 y/l - i 2 + fi, f" 1 = 2 - n + 2 Vl-t 2 -/^ . (28) 

If 1/2 < £ < 1 then ^ and £2 are complex conjugate of each other, whereas if0<£<l/2 
they are distinct real numbers. Also notice that 

6+er 1 = 4 + 2/i, £ 2 + £ 1 = 4-2 f jL. (29) 

Indeed, in order to verify formula (J25|) it suffices to remark that the right hand side equals 

(£1 + -V- y- l m + & - y - I/" 1 ) = 16 - V - 8(y + y- 1 ) + (y + y- 1 ) 2 

which is the same as the left hand side. 

Proposition 4.2 For < t < 1, we have 

E{<j>) = (l-eija-eDa-e^a-^oa-^^-V), (30) 

G(*-V) = 77^- (31) 

Proof. First of all notice that, by Lemma I4~TI = e ai/2+bQ and cr -1 = e a2/2+bQ with 

ai = — - log(l — 2t cosx + t 2 ) + iri, a 2 = ai — log a. 
We apply Corollary 12.41 twice (with bQ instead of Q) and obtain 

E{4>) = exp (2trace(#(ai)iJ(ai) - #(a 2 )#(a 2 )) Wcr~V)- 



Introducing a± = log(l — 2tcosx + t 2 ) and a 2 = log(t 2 + sin 2 a; + sin 4 x), we have a\ 
—a.i/2 + 7t~i, a 2 = a\ + ai + a 2 /2 = [oti + ct 2 )/2 + ni (see ()24|) ). and thus 

exp ^2 trace(i?(ai)if(ai) — H(a 2 )H(a 2 )) 

= exp ( - ~ trace(iy(ai)if(a 2 ) + H(a 2 )H(a 1 ) + H(a 2 )H(a 2 ))^ . 
With the substitution 2 = e lx we can write 

= log((l - tz){l - tz' 1 )), a 2 (z) = log((z 4 - 8z 2 + (14 + 16t 2 ) - 8z~ 2 + z' 4 )/l6). 
Since we can decompose 
z A - 8z 2 + (14 + 16t 2 ) - 8z- 2 + z- A = ^ 2 \1 - ^z 2 ){l - i 2 z 2 )(l - £iz~ 2 )(l - ^z~ 2 ), 
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we have the Fourier coefficients (k ^ 0) 



[«l]fe = i a ^hk = 1 2 , [a^Wl = 0. 

We obtain 

tracetf(ai)#(a 2 ) = traced (a 2 )H(a 1 ) = V t 2fc - ^ 2 = - log((l - t 2 £i)(l - t 2 ^)) 

k=l 

and 



tracetf (a 2 )tf (a 2 ) = £ ^** 2j = -2 log((l - - £ 2 2 )(1 - 66) 
fe=i 

This implies 

E(d>) = (1 - £ 2 )(1 - £ 2 2 )(1 - 66) 2 (1 - t%)(l - t%)E{a- l <t>). 
In order to compute G(cr _1 0) notice that 

det(cr~V) = ^ ^ = (t 2 + sin 2 x + sin 4 x)(l — 2t cosx + t 2 ) = exp(a 2 + cm)- 
a 2 

The arithmetic mean of oti is zero while the arithmetic mean of a 2 is log(l/(1666)) as can 
be seen from the factorization. □ 

Clearly, the previous proposition reduces the computation of E((f)) to the computation 
of E(ip) where 

/ p q\ 

i> := a~ l cj> = . (32) 

\q p J 

Since the Fourier coefficients ipk vanish if \k\ > 3, the constant E(ip) can be computed from 
Theorem 12 . 71 with n = 3: 

Ety) = G(^) 3 detT 3 (^- 1 ) . (33) 

We know in addition G{ip) = G(a~ l 4>) so all that remains is the computation of the determi- 
nant of a block 3x3 Toeplitz matrix whose symbol is the inverse of ip -1 . The next section 
is devoted to this task because even the determinant of a 3 x 3 block matrix (which is in 
fact 6 x 6 in size) can be difficult to compute. 
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5 The computation of detT 3 (^ x ) 



The symbol ip 1 is given by 



ip 1 = i] 





(34) 



where rj is the even function 

rj = (pp-qq)' 1 = (det(cr~V))^ 1 = (1 - 2t cosx + t 2 ) _1 (t 2 + sin 2 x + sin 4 a;)" 1 . 

Using the fact that b is an odd function, it follows that the block Toeplitz matrix T 3 (?/>~ 1 ) 
has the following structure: 

\ 



a 





ai 


61 


a 2 


&2 





a 


-61 


a_i 


-62 


a_ 2 


a_i 


-h 


a 





ai 


61 


h 


ai 





a 


-61 


a_i 


a_ 2 


-62 


a_i 


-61 


a 





b 2 


a 2 


61 







a 



V 

What we will next attempt is to show is that the special form of this matrix allows us to 
reduce its determinant to that of a 2 x 2 scalar matrix. The first thing to notice is that if 
Ci = det A and c 2 = det B then the above matrix has the form 

a I 2 A B 
C\A~ X aoh A 
c 2 B~ l ciA^ 1 a I 2 

Now if we factor out the a terms and multiply the top block row from the left by CiOq 1 A~ 1 
and substract this from the second block row we are left with 

h a x A a^B 

(1 - cia 2 )/ 2 clq X A - da^A^B 
c 2 a 1 B- 1 Cl a^A- 1 I 2 

Except for the factor a® this has the same determinant as the one of interest. Multiplying 
the top block row on the left by c 2 a 1 B^ 1 and then subtracting yields in the same manner 
a determinant equal to that of the matrix 



a, 1 A 



c±a 



(1 - cia 2 )/ 2 
i A -i 



c 2 a^B- x A 



a Q l B 
i^A - ciOq 2 A~ 1 B 
(1 - c 2 a 2 )I 2 
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The determinant of the above is given by 

det((l - ci% 2 )(l - c 2 a 2 )I 2 - (ci% 1 A~ 1 - c 2 % 2 B~ 1 A)(a 1 A - aa^ 2 A~ 1 B)) 

or 

det((l - 2cia 2 - c 2 a 2 )I 2 + c 2 a 3 B^A 2 + c\a z A' 2 B) . 

Now any 2x2 matrix of the form (det M)M~ l + M is a constant times the identity. This 
is easy to see if one thinks of the inverse of the matrix. In fact if the matrix is given by 

M = 

then the constant is mi + m 2 . The matrix 

c 2 B- x A 2 + c\A- 2 B 

is of this form and it follows that the matrix given in the last determinant expression is a 
constant times the identity. Recalling the definition of the terms A and B we see that this 
constant is 

a o 3 ( a o ~~ 2cia - c 2 a + a_ 2 (a 2 - b\) + 2b x b 2 {a x + a_i) + a 2 (a 2 x - b\)). 

Since we originally had a factor of the 3 term cancels (since this comes from a 2 x 2 
matrix). Thus, at this point, we have shown that 

detT 3 (^ _1 ) = A 2 (35) 

with 

A := al - 2c 1 a - c 2 a + a_ 2 (a 2 - b\) + 2bib 2 (ai + a_i) + a 2 (a 2 _i - b\). (36) 



Now we will investigate the Fourier coefficients. Going back to ()23|) and the functions 
a and b are given by 



He") = , 2+ .f' . 0(0 = E^' + t''' . 4 (37) 

r + sin i + sin x (t — e " ) (i^ + sm i + sm x) 



With the substitution z = e tx the function 6 can be written as 

b( Z ) : 



1 — z) 



z 4 - 8z 2 + (14 + 16t 2 ) - 8z- 2 + z~ 4 

8^2 (a" 1 ~ z) 

(1 - ^ 2 )(1 - 6^ 2 )(1 - iiz~ 2 ){l - &- 2 ) 
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The Fourier coefficients of this can be easily computed and a little algebra shows that 
b 2 = and 

1 (i + £0(1 + 6X1 -66)' 

With 62 = we also can simplify the constant (J3l)j) to 

A = <2q — 2(ai<2_i + b\)dQ — 020.-2^0 + 0-2(^1 — b\) + 02 (ai^ — & 2 ). (38) 
The function a can be rewritten as 

8t(z + z- 1 ) - A(z + z- 1 ) 2 



a(z) 



(t - z- 1 )^ 4 - 8z 2 + (14 + 16t 2 ) - 8z- 2 + z~ 4 ) ' 



Using the factorization employed for finding the coefficients for b it can be shown in a messy, 
yet elementary computation that the following holds: Let 

° := (i -&&)(£, -6) (39) 

and 

x 2 -4x-l 1 / 2 2 \ 

(x) : ~ (1 - t 2 x){l - x 2 ) ~ T^¥x~ [TTx' ~ T^x~ ~ ) ' { ' 

(1 - 2t 2 )x 2 + (-2 - 2t 2 )x + 1 . 2 

/(x) := z s-V; ^ = fc(a:) H . (41) 

V ; (1 -t 2 x)(l -x 2 ) y 1 1 -x y ! 

Then the Fourier coefficients of a are given by 

a = ta(^fc(^i) - 6&(6)), 
fll = a(Z(£i) -/(&)), 

a 2 = ta(k(^) - k(g 2 )), 
With this notation our constant (|3*H1) becomes 



A = a 3 (^ - 6) 2 (^i*(Ci) 2 *(&) - t 3 6M6) 2 M6) 
+tk{t 1 )ifa)* - ^(6)^i) 2 + 4(1 + £i)- 2 (i + 6)- 2 (^(6)(6 + 1) 2 - f*(6)(6, + 1) : 

Simplifying this expression is the last remaining task. 

For simplicity we will assume that 1/2 < t < 1. In this case the roots £ and £2 ar e 
conjugate and \i is complex. The case < t < 1/2 is similar, but we omit the details. 
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The answer however is valid for all < t < 1 because of analyticity of the expressions with 
respect to t. We begin by proving that the expression 

t^ik^fm) - t 3 6M6) 2 M6) + t*(&)l(6) a - f*(6)l(6) a 

4*M6) 4tfc(6) ^ 42 ) 
+ (i + 6) 2 (1 + 6) 2 

is equal to 



8/i 



(43) 



(l + £i)(l + 6)V^ 
where 

w:=(l-t 2 6)(l-t 2 6)>0 (44) 

First recall the identities 

ei + er 1 = 4 + 2/i, 6 + ^ = 4-2/1, /i 2 = 1 _ 4t 2 (45) 

from which it is easy to conclude that 

(6 - IKS" 1 - - 1)^ 1 - i) = let 2 . (46) 

In order to evaluate the expression in (J42)) we first simplify 

P£Mi) 3 kfo) - t 3 6M6) 2 M6) + *M6)/(6) 2 - t*(6)J(&) a - (47) 

Using ({1Q}) and (jUj) the expression 
becomes 



^*(£0 + ^—H = *6M6r-^(6) 



^ vsi/ ^ vsi/ 1-6 (i-6) 2, 

A little more algebra shows that this is the same as 

, n / 2 2 \ 4t 

-rrr - 1 — r + 1 



1 + 6 1-6 y (i-e 



V2' 



Since we have only used the identity for in terms of the above is also true if we 
replace £1 with £2- Using this last above equation (for both £1 and £2 terms) we have that 
(|4Tj) is given by 

*fc)*(6) (-T^T -T^ + T^ + T^)- + ^ 



1 + 6 1-6 1 + 6 1-67 (1-6) 2 (1-6) 2 
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which is of course 

itk( nk(n( \ 4 ^(6) 4^(6) 

We now proceed to simplify the above. From (|45j) it is easy to show that 

(1 - M 6) 2 = 4^(1 - t 2 6), (1 + /i6) 2 = 46(1 - t 2 6), 
which using ()46p implies 



= (i-^i) 2 (i + ^ 2 ) 2 = (i-^i) 2 (i + ^ 2 ) 2 (6-i) 2 fe-i) 2 Uq) 
1666 ~ (16660 2 ' 

Moreover, using 2(1 — /i6) = 1 — £ 2 + 46 and 2(1 + /i6) = 1 — £ 2 + 46 (which also follows 
from (jEJ)) we conclude 

M6) = M6) = (i + ^ 2 )(i-6 2 )- 

With this we have that 

mm2) ~ (i-^i)(i + ^ 2 )(i-6 2 )(i-a)' 

Using (J4UJ) it then follows that 

4 

tv^7(l + 6)(1 + 6) 

Hence the term (|4"7j) (which is also (|48)l) is 

16(6 2 -6 2 ) 4tA;(6) + 4tA:(6) 



v ^(i + 6)(i + 6)(i-6 2 )(i-^ 2 ) (i-6) 2 (i-6) 2 ' 

and the expression for the constant (|4"2j) is the above plus the additional terms of 

4tfc(6) 4tfc(6) 
(1 + 6) 2 (i + 6) 2. 

In other words we still need to evaluate 

i6(6 2 -£ 2 ) 



v ^(i + 6)(i + 6)(i-6 2 )(i-e 2 2 ) 

4tfc(6) 4tfc(6) 4tfc(6) 4tfc(6) 



(1-6) 2 (i-6) 2 (i + 6) 2 (i + 6) 2 ' 
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(50) 



If we combine the terms 

4tfc(£0 4tfc(&) 



(l-£ 2 ) 2 + 

we have 

8tfc(£Q(l + £ 2 2 ) 

(i-£ 2 ) 2 (i + £ 2 ) 2 

which becomes after using the identity for fe(£0 

-64^(1 + £ 2 2 ) 



(i-^i)(i-ex 2 )(i-6) 2 (i + 6) 2 ' 

Using (J4T?|) and (|43|) this becomes 

-4(^ X + + ^2) "8(2 -//)(! + ^£2) 



+ £0(1 - 6)(i + 6) 2 v^(i + + - £ 2 2 ) 

or 

-16 + 4//(2 - 46 + 4//£ 2 ) 



^(i+eoa+e^a-e: 



and finally 

-16 4u 
+ 



v^J(i + Ci)(i + 6)(i - £f) v^(i + £0(1 + 6) ' 

We still have another term same as the above but with the £1 and £2 switched, \x replaced by 
—\l and an additional sign change. However when we add this new term to the above and 
then substitute in (|5U|) we have that our constant is 



v ^(i+e 1 )(i+6)(i-ei 2 )(i-a) 

-16 16 

■ 

8/2 



^(1 + 6)(i + £ 2 )(i - e 2 ) v^(i + eoci + - 



(51) 



^(l + £i)(l + £ 2 ) 

which nicely reduces to 



v^(i + 6)(i + 6) 

as promised. Thus the constant (|55|) is 

8^(6 -£ 2 ) 2 466 r , ^ 

A = —=-. — with a = — -. 52 

v^(i + £0(1 + 6) (l-£i£ 2 )(£i-£ 2 ) 1 ; 

Now we can put the final pieces together. 
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Theorem 5.1 For < t < 1, we have 



2t(2 + t 2 ) + (l + 2t 2 )v / 2n r 
Proof. We have from Proposition 3.2 and formulas (jH2J), (|33|). and (JH2J) that 

£(</>) = (l - - 6 2 )(i - 66) 2 (i - - * 2 6)(1666)- 3 

8/ ^ 3 (6-6) 2 x2 
^1 + ^(1+6) 



This can be reduced to 

4 5 tV^ 2 4 



(i-6 2 )(i-6 2 )(i-66) 4 (6-6 



This uses only elementary algebra and formula 

To evaluate the remaining terms we note (J45|) implies that 

and 

(1 + £0(1 + ^(l + 6)(1 + G" 1 ) = (6 - 2/i)(2 + 2/i) = 16(2 + t 2 ). 
From the last equation together with (}4"Tj) we can conclude that 



(1 - - £*) = 16^(2 + ^16. 
Thus from the line above and (fo"4"j) it follows that 

EM = . 

To complete the computation notice again ([45)1 . hence 

(6 + 6 _1 )(6 + = (4 + 2/x) (4 - 2//) = 12 + 16t 2 

Also we know 

(ST 1 - 6X6" 1 - 6) = 16^(2 + t 2 ), 

so we can check that 

66 + 1/66 = 6 + 8t 2 + 8V* 2 (2 + t 2 ). 

This yields that 

66 i 



(1-66) 2 4 + 8t 2 + 8v/t 2 (2 + i 2 ) 
From this the theorem follows. 
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6 The asymptotics for more general parameters 



In the previous sections we have shown that 

det M n ~ E{4>), n -> oo, (55) 

if < t < 1 and we have identified the constant E{$). The goal of this section is to show 
that (jB3j) holds for all values of parameters t belonging to the set 



n 



[teC : Re (t) > J 



under the assumption that the entries of M n are continued by analyticity in t onto Q. 
Indeed, if we recall the definition of the entries of M n in terms of Rk and Qk (as stated in 
the introduction) one only needs to observe that the denominator of the fraction under the 
integrals is nonzero as long as t 2 is not real non-positive number, i.e., if Re (t) ^ 0. 

Since we have related M n to T n (</>) by simple row and column operations (see (J2J) or (|17|)) 
we are going to elaborate on T n (<j)), the entries of which also depend analytically on t. Recall 
that is given by 

c t d t 



with 



ct(0 



d t c t 



t cos x + sin 2 x 



d t (e lx ) 

Moreover, we introduce the function 

et(e lx ) = 



' e -ix _ t) Vt 2 + sin 2 x + sin 4 x 
sin a; 



Vt 2 + sin 2 x + sin 4 x 



t cos x + sin 2 x 



(e~ ix - t)Vt 2 + sin 2 x + sin 4 x e ix - t' 

Proposition 6.1 For each it G the functions d t and ef belong to B, and the dependence 
on t is analytic. 



Proof. The statement concerning the function t is east to see. Indeed, the function 
sin 2 x + sin 4 x belongs to B and has its spectrum equal to [0, 2]. By functional calculus we 
can define the reciprocal of the square-root of t 2 + sin 2 x + sin 4 x whenever t E C\ [— i \/2, i\/i] 
and also at t — oo. This domain (extended with the point at infinity) is simply connected 
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and we can choose an analytic branch. Cleary, we choose the branch in such a way that the 
square-root is positive for real, positive t. 

In regard to the function e t the assertion is obvious except for \t\ = 1. We rewrite this 
function as 

, (t cos x + sin 2 x) — Vt 2 + sin 2 x + sin 4 x 

e t (e ) = - . 

(e~ ix - t)Vt 2 + sin 2 x + sin 4 x 

We are going to show that term e~ lx — t cancels with a term in the numerator for \t\ = 1, 
Re (t) > 0, and that the resulting function belongs to B and depends analytically on t. To 
see this rewrite the numerator (with the substitution z = e %x ) as a quarter times 

n(t, z) = 2t(z + z' 1 ) -(z- z- l f - Vz 4 - 8z 2 + 14 + 16t 2 - 8z~ 2 + z 4 . 

Allow, for a moment, z to take complex values in a neighborhood of the unit circle. More 
specifically, it is easy to see that for each e > there exists a 5 > such that the function 
under the square-root is nonzero (and of course analytic in both z and t) whenever Re (t) > e, 
\t\ < 1 + e, 1 — S < \z\ < 1 + 5. Denote the set of all (t, z) satisfying these conditions by U £ . 
Hence n(t, z) is analytic on U £ . Consider the subset M £ — {(t, z) G U e : z — 1/t}. Then 

n(t, r 1 ) = {t 2 - r 2 + 4) - ^(t 2 - t~ 2 + 4) 2 . 

For t = 1, keeping in mind the proper choice of the branch of the square-root, this equals 
zero. Hence by analytic continuation and since M e is connected and the function under the 
square-root does not vanish on M e it follows that n(t,z) = for all (t,z) G M e . Hence 
n(t, z) = (z" 1 — t)m(t, z) with some function m(t, z) which is analytic on U e . From this the 
assertion follows. □ 

Since the function d t depends analytic on t G Q, the entries of the matrix T n (d t ) are also 
analytic in t. In regard to the matrix T n (ct) we decompose for < t < 1 

T n (c t )=T n (ef) + K+ t , (56) 

where K^ t is the n x n Toeplitz matrix with entries 




if j < k 
if j > k. 



This decomposition holds since K n>t is the Toeplitz matrix with the generating function equal 
to (e~~ lx — t)^ 1 for \t\ < 1. The right hand side of (jo^j) is well-defined and analytic for t G tt, 
and hence we can define 

B ntt :=T n (et)+K+ t , 
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which is the analytic continuation of T n (c t ) onto Q. From this it is easy to see that the 
analytic continuation of T n ((f>) onto t G fl is the matrix 



B n ,t = 



B n ,t T n {d t ) 



Tn(d t ) T Bl, 



where A T stands for the transpose of A. 
Theorem 6.2 For each t G fi 7 we have 

lim det B n t = E t with E t = (57) 

Moreover, the convergence is locally uniform on Q. 

Proof. We make the decomposition = 0+ + O with 

et d t \ ( (e'^-t)- 1 

0+= | , 0o(O = 

d t et J V (e^-t)- 1 



In view of the definition of B Utt we obtain 

B n , t = T n (<f> + ) + K:, K+:= 

Now let 



Kt 



o (Kt) 



T 



1 - te lx 
1-te 



0+(O = 

A straightforward computation shows that 

T n {Q + )Kt = T n (0+0 O ), (©+0o)(O 





e~ ix 



Moreover, 

T„(0+)T n (0+) = T n (0+0+) - P n tf (0+)tf (0+)P„ - W n H{Q + )H{<t> + )W n . 

Adding the above two equations, taking the determinant and noting that detT n (0 + ) = 1, it 
follows that 

det B n , t = det (V n (0) + P n KP n + W n LW n ^j , : = 0+0+ + 0+0 o = 0+0. (58) 
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The function belongs to J3 2x2 and depends analytically on t G Q. Moreover, det = 1, 
hence G QB 2x2 for each t E Vt. For < t < 1, Lemma E~T1 implies that G £i£> 2x2 , 
and since Q is connected this holds for each t E Q. The operators K and L are trace class 
operators (in fact, rank one operators), which also depend analytically on t E Q. Hence all 
the assumptions of Theorem 13. II are fulfilled, and applying it yields the desired assertion. □ 



Appendix A: The basic identity 

As mentioned in the introduction, the block matrix considered by the authors in [Hj is of the 
form 

Tl Q 

M n = | J (59) 

Q n 

where TZ andQ are n x n matrices with entries given by 

TZ jk = 2(-l f k -^R k ^ +1 + 6(j - k^-*- 1 
Q jk = 2i{-lf^^Q n+1 _ 3 _ k 

(1 < j,k <n), and the expressions R k , Q k , and 6{k) are defined by the following. 
For even k 

i r r cosy cos (&£ + ?/) 

k = sTT / / 2 1 2 7T2 2? —-AV dx i 

^tx 1 J_ w J_ w cos z x + cos z y + t z cos z {x + y) 



and Q k = 0. 
For odd k 

and 



Y r-K COS X COS kx 

Qk = -z-j / / — — — 2 . , 2 27 — : — zdydx, 

8ir z J_ n J_ T cos z x + cos z y + V cos z {x + y) 

t H f w cos(x + y) cosikx + y) , 

Rk = FT / / — 2~^ 2 _, - 2 i ir —^dydx. 

6ir z J_ n J_ % cos z x + cos^ y + t z cos z {x + y) 

The expression 8(k) equals 1 for k > and otherwise. 

We begin with the following lemma which describes the above coefficients as Fourier 
coefficients of certain functions. This is our starting place, since once this is done, we use some 
simple matrix algebra to express the determinant of the above matrix as the determinant of 
a block Toeplitz matrix. Throughout what follows we will assume that < t < 1 because 
the second part of the matrix 1Z becomes unbounded in the limit n — > oo in the case t — 1. 
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Lemma 6.3 Consider the functions 

1 ' Air J_ n cos 2 (x-n/2) + cos 2 y + t 2 cos 2 (x + y-ir/2) U ' 

T( e lx ) = - / V dy. 

An J_ n cos 2 (x — 7r/2) + cos 2 y + t 2 cos 2 (x + y — it/2) 

These functions have Fourier coefficients, S k and T k , satisfying Sk + T k = (— l)^~^R_ k+1 
where R k is as defined as above. 

Proof. Replacing x by x + it in the integral defining S shows that S(e l( - x+n ^) = S(e tx ) and 
hence that S k = for odd k. For k even we have 



/7T PIT 
-IT J —IT 



cos(x + y- 7T /2)e i(x+ y- 7T We- ikx 
87T 2 J_ n J_ w cos 2 (x — 7r/2) + cos 2 y + t 2 cos 2 (x + y — 7r/2) 



dydx. 



If we replace x with x + n/2 then this integral becomes 

t(_i)$ r r cos^ + ^e^-^- 1 ^) 



/7T /»7T 
^ J —-jr cos 2 x + cos 2 y + t 2 cos 2 (x + y) ^ X 

k 

t(—\p f n f n cos(x + y) cos(y — (k — l)x) 

cty ctoc . 



87r 2 J_ 7T J_ n cos 2 x + cos 2 y + t 2 cos 2 (a; + y) 

The last equality follows since the inner integral as a function of x is even and this is exactly 
for even fc. If we replace x by x + 7r in the definition of T shows that the 
even coefficients vanish. A very similar computation to the above shows that for odd k, 
T k = (— l)' _ 2'i?_fe+i and thus combining the two cases we have S k + T k = (—l)^^R_ k+1 for 
all k. □ 

Lemma 6.4 ForO < t < 1, the matrixlZ is annxn Toeplitz matrix with symbol 2(S+T)+U 
where S and T are defined in Lemma 1 and U(e lx ) = e -L_ t - 

Proof. This follows directly from Lemma 1 and the fact that the Fourier coefficients of 
-h- t aret fc -^(A;). □ 

Lemma 6.5 Define 

Vie-) = (_l)[fl+iJ- f C ° S ^ ~ *M dy. 

v 1 v 1 Air J^ n cos 2 (x -ir/2) + cos 2 y + t 2 cos 2 (x + y - ir/2) 

Then V k is zero for k even and for k odd we have V k = i(— 1) t" + 2 k ^Q k 
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Proof. To see that Vk is zero for k even, replace x by x + tt as before. For the remaining 
equation we have 

k 1 J 87r 2 J /_ 7r> /^cos 2 (a:-7r/2) + cos 2 i/ + t 2 cos 2 (a; + i/-7r/2) ^ 

Making the substitution x + 7r/2 for x and simplifying yields the result. □ 

Our next step is to take the functions given by these integal representations and simplify 
them. To this end, notice that the denominator in all three integrals, 

cos 2 (x - it/2) + cos 2 y + t 2 cos 2 (re + y - n/2) 

can be written as 

sin 2 x + cos 2 y + t 2 sin 2 (a; + y) 

which is in turn 

-({2 + t 2 -cos2x) + (1 -t 2 cos2x) cos2y + t 2 sin2a;sin2|/). 

If we combine the integrals for S and T and use the obvious trigonometric identities we see 
that the numerator is of the form 

-{{t - e lx ) - {te 2lx + e ix ) cos2y - i(te 2ix + e ix ) sin2y) 

Thus combining S and T we see that we have an integral of the form 



2?r 



D + E cos 2y + F sin 2y 
An J A + B cos 2y + C sin 2y 



dy, 



where the coefficients depend on x and t and are given by the above two equations for the 
numerator and denominator. 

It is fairly straightforward to evaluate such an integral. First we replace 2y by y which 
results in no change because of periodicity. Second we define 

B C 
cos0= =, sm< 



V5 2 + C 2 ' VWTc 2 

and change variables using y = z + 0. At this point we have an integral of the form 

1 f 27T D + G cos z + H sin z , 

-dz. 



An In A + K cos 2; 
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The contribution from the sin-term is zero and the rest can be easily evaluated since it 
reduces to one of the form 

dz 2tt 



o 1 + 0. COS Z y/l - a 2 

providing that \a\ < 1. Using this equation we can compute the integral for S+T. In the same 
way we are also able to evaluate the integral for the function V. This computation is even 
more straightforward and we leave the details to the reader. Notice that the denominator in 
the defining integral is the same as for S and T and the numerator only depends on x. 

Lemma 6.6 We have 

t cos x + sin 2 x 1 

o -\- 1 
V 



2(t - e~ ix )Vt 2 + sin 2 x + sin 4 x 2(t 

(_l)K2]+l sinx 



+ sin 2 x + sin 4 x 

Hence the generating function for the Toeplitz matrix 1Z is given by 

t cos x + sin 2 x 



2(S + T) + U 



+ sin 2 x + sin 4 x 



Notice finally that the (J, fc)-entry (1 < j,k < n) of the Q matrix is given by twice the 
n + 1 — j — k Fourier coefficient of V. 

Now that we have identified the entries in each of the blocks of the matrix (jf>9j) as Fourier 
coefficients of certain functions, it is fairly easy to see how the determinant of (|59jl can be 
computed as the determinant of a block Toeplitz matrix. As a preliminary step, we define 
an operator W n on the finite-dimensional complex vector space of dimension n by 

W n (at, ...,a n ) = (a n , ...,a 1 ) 

and denote the identity operator on the same space by I n . Multiplication of a matrix by W n 
on the right results in changing the (j, fc)-entry into the (j, n + 1 — fc)-entry and on the left 
changes the entry to the (n + 1 — j, fc)-entry. Also note that det W% = 1. Hence the block 
matrix product 

I n \ / K Q \ / I n 




o w n J \ q n J \ o w n 

has the same determinant as the one given in formula ()59|) and thus we can replace it with 
the product. The advantage to this is that since 

In0\/KQ\/I n 0\ ( K QW n 




W n ) \ Q K j \ W n ) \ W n Q w n nw n 
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the (j, fc)-entry of QW n is now the k — j Fourier coefficient of the function V and the (j, k)- 
entry of W n Q is j — k and thus the upper-right and lower-left matrices are now Toeplitz in 
structure. The matrix W n VW n still remains Toeplitz, but with symbol (2(S + T) + U)(e lx ) 
replaced by (2(S + T) + U)(e~ ix ). 

The end result of this is that we have shown that (for < t < 1) 



detM n = detT n (0) with = 




where 

t cos x + sin 2 x 

[e~ %x — t)Vt 2 + sin 2 x + sin 4 x 
sinx 

+ sin 2 x + sin 4 x 

c(e lx ) = c(e~ tx ), d(e lx ) = d(e~ tx ) Notice that dropping the sign when passing from V to d 
does not influence the value of the determinant. 
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